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Abstract 

Many methods for reducing and simplifying differential equations are known. They provide 
various generalizations of the original symmetry approach of Sophus Lie. Plenty of relations 
between them have been noticed and in this note a unifying approach will be discussed. 
It is rather close to the differential constraint method, but we make this rigorous basing on 
recent advances in compatibility theory of non-linear overdetermined systems and homological 
methods for PDEs. Q 

Introduction to reduction methods 

Exact solutions play an important role in investigation of equations of mathematical physics. 
Even though usualljo it's only a small portion of solutions one can find explicitly, they represent 
significant physical phenomena and are often applied in practice. The first explicit methods have 
origins in the works of Euler, Lagrange, Cauchy, Hamilton, Jacobi and Monge. But the most 
essential contribution was made by Sophus Lie, whose continuous transformation groups changed 
the landscape of the theory of differential equations. These transformations are called nowadays 
symmetric^ and can be computed with the help of most popular symbolic software packages. 

If the symmetry algebra is known, is sufficiently large and has an appropriate Lie structure, 
then the given differential equation (or system) can be integrated completely. This works both for 
ODEs and PDEs. But usually the algebra is not bi^, so other methods should apply. 

There is a number of classical methods, like intermediate integrals, differential and nonlocal 
substitutions (e.g. Laplace, Backlund and Darboux transformations) and conservation laws. In 
more recent time a lot of other methods for obtaining exact solutions were invented. The whole list 
is large, but we name a few: higher and inner (=internal) symmetries [KLV] . recursion operators, 
nonlocal and ghost symmetries, nonclassical symmetries and direct reduction, inverse scattering, 
bi- linear Hirotha method, Sato theory and so on, see e.g. [BKi ICRCl ISCLj . 

Many interrelations have been made explicit, but the unifying theory does not exist so far. In 
this paper I will try to present a moderately general viewpoint for all local theories (nonlocal ones 
follow the same scheme, but additionally need a covering, i.e. another system of DE implying the 
given one as a differential corollary). I will also try to explain why it is cumbersome to encompass 
all features into the framework of one theory leaving it reasonably efficient. 



0.1. Differential constraints. 

This method is pretty old (it dates back at least to the time of Lagrange) and is quite obvious: 
given a system £ — {Fi ~ ■ ■ ■ — Fr — Q\ one overdetermines it with other differential equations 
Q — {Gi = • • • = Gs = 0} and investigates the joint system £ r\Q = {Fi = 0, Gj = 0}. Clearly 

iMSC numbers: 34A05, 35N10; 58A20, 35A30. 
Keywords: Compatibility, differential constraint, symmetry, reduction, multi-bracket, solvability. 
^Except for the important class of completely integrable equations. 

^They are distinguished to be point or contact symmetries, though the latter received less circulation. 
''For linear system it is huge, but the bulk consists of shifts by solutions and the quotient is usually moderately 
small. 
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Sol{£ O Q) C Sol{£) (here and in what follows by Sol(5) we understand local or even formal 
solutions, though global will work as well). 

Remark 1. This procedure is the direct analog of the standard tool in algebraic geometry: In order 
to investigate an algebraic variety X C it is bi-sected by another variety Y and the totality of 
intersections X (lY uncover X ( most often Y is taken in the simplest way, like an affine plane 
pfc pAf^ sometimes more complicated variety is needed). In the case of differential equations 
the section Sol{Q) may be even more complicated than the original Sol(£), but the task is to make 
Sol(£:n5) tractable. 

We will suppose that the original system £ is compatible (=formally integrable). Let us call 
the constraint Q compatible differential constraint if the joint system £ n Q is compatible as well 
(NB: Q is not required to be compatible!). 

This method was introduced into practice by Yanenko jY] and also as a generalization of S.Lie 
approach by Ovsyannikov [Ov| . see also Bluman and Cole |BC| . though relations between these 
were unexplored for some time |Me| . Re-invented by the name of side conditions in [ORj . the 
method was much popularized later by Olver |01j . In fact, it was in this latter papers that the 
most cited statement was made (we change it a bit and state the claim in the ultimate form): 

All reduction methods are partial cases of the method of compatible differential constraints. 

In particular it was shown in ^Olj that many known methods such as classical and non-classical 
symmetry reductions, direct method of Clarkson-Kruskal, Galaktionov's separation method and 
some others are covered by the compatible constraints. 

This is an important statement and some remarks should be of order. 

1. The result is demonstrated for a scalar single second order PDE on the plane only and even 
though it was claimed that "the results can be straightforwardly generalized to arbitrary systems" , 
it is neither clear to what degree (see 2 below which shows one must be precise) nor how (direct 
generalization of this proof will lead to the trouble even with a single scalar PDE of order 3 on 
M'^; consult [KLiIIKLTI for a similar but rigorous statements to appreciate which technique should 
be used for quite general, though still not arbitrary systems). 

2. The statements in |0R1 lOlj are rather vague, but without specific genericity conditions the 
results turn wrong. For instance the following is a combination of Proposition 1 and Theorem 2 
from [01] : A first order differential constraint G = reduces the second order PDE F = to an 
ODE iff the system F — 0,G = is compatible. The trivial counter-example to this is F = H oG, 
which is compatible with G for any differential operator H (if the operators arc non-linear we 
should rather use linearization in the composition: F — in o G). 

3. The situation becomes worse in higher order. For instance Darboux semi-integrable systems 
can be treated as compatible pairs F = 0, G = of scalar PDEs on the plane with 1 common 
characteristic, which surely declines possibility of direct reduction to an ODE. For higher order 
systems on R" the situation becomes completely untractable if no restrictions are imposed (the 
whole algebraic singularity theory imbeds into this problem!). 

Lacking a general proof and due to these remarks we regard the above claim as a conjecture, 
whose precise form should be specified. The proof will be given in the subsequent sections. 

0.2. General statement. 

Reductions of systems of PDEs to systems of ODEs lead to final dimensional families of solu- 
tions. Actually this can be taken as a definition of reduction: Any stratified manifold S C Sol(f ) 
will be called a reduction to ODEs. This in particular covers the case when ODEs have singular 
solutions because these should be considered as components of S of different dimensions. The 
claim formalizing the one of the previous section is now the following: 

Theorem 1. Any reduction S C Sol(£) can be obtained via the method of compatible differential 
constraint. 
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This statement is not too rigorous either. The main problem is the space Sol(£) that is not 
defined due to the lack of existence and uniqueness theorems in PDE theory. Some solutions to 
this exist. One is to consider analytic systems and analytic solutions or to boil everything down 
to the formal theory (this requires a lot of formalism but is possible). Another one is to restrict to 
classes of smooth PDEs, which have existence and uniqueness theorems (together with regularity 
theorems so that a-posteriori we do not need to leave the class of smooth functions). We would 
prefer however to consider the stratified manifold S abstractly and claim that all its points are 
(loc.) solutions of £. There are standard ways to impose the topology on it. 

Proof. We can fix a setup so that S C C°°(7r) for a bundle tt : E ^ M. It is possible to do 
globally due to Whitney theorem for stratified manifolds, but local representation will suflice for 
us because we'll be patching them later on. 

Consider all (non-linear) differential operators G G diff(7r, 1) such that G[u] — for any m e 5. 
This is a set Y{S) D Y{£) = {Fi, . . . Fr). The space Y{S) is a differentially closed algebra. Indeed 
if Gi G Y{S) vanish on section u G 5, so do their linear combinations and left-compositions with 
other differential operators. 

Now we claim that Y{S) is finitely generated over diff(l, 1). In formal theory this is so due 
to Noetherian property (one can consult IKL5I for the framework of formal theory for non-linear 
operators, when modules should be considered over the algebra of "^^-differential operators). 

In analytic context this is generally wrong, but we can assume that Y{S) is a prime ideal and 
then this becomes the fundamental theorem of Ritt-Kolchin theory [Ril IKoj (we refer to [BJ IMaj 
for the theory of D-modules and applications to our context). 

In smooth case one should act as follows: a neighborhood of a stratum in S is parametrized 
by constants Ci, . . . , q. Taking the space of all parametrized derivatives D„u{ci, . . . , q) the im- 
plicit function theorem makes it possible to specify the defining functions of the d-dimensional 
submanifold with parameters This finite set can always be specified to begin with Fi, so it is 
Fi, . . . , Fr, Gi, . . . , Gs- The second part of this collection defines Q. 

Note however that in any of the above approaches it is not possible to distinguish the whole 
manifold S at once. Overdeterminations Q play the role of local charts and as in differential 
geometry we need a finite or countable collection of such. This is especially clear with application 
of implicit function theorem, but is necessary for other cases as well. Finally we patch the local 
charts and so S is defined by the compatible differential constraint Q. □ 

Even though on the level of ideas this proof is clear, the reduction methods have quite precise 
implementation techniques and algorithms. A prompt objection to the above is that it is not very 
constructive. But this will be compensated in the next section, where we explain how standard 
methods imbeds into compatibility theory and make specifications. 

Remark 2. The above reduction to ODE, in light of the analogy from Remark]^ is similar to 
the section of X by transversal varieties Y of complimentary dimension, codiml" — divaX, when 
dim(A' C]Y) — Q. Note however that with the above approach we always obtain not only solutions 
S , but more generally its D-Zarissky closure S . 

Higher- dimensional strata of So\{£) are parametrized by a functions of p variables (a is the 
functional rank and p the functional dimension o/Sol(£), we refer to IKL^^ for the discussion). 
Thus we can also make a more general reduction of p or keep it and reduce a. 

Note however that several small reductions are not equivalent to one big, because several steps 
can be impossible, see 'C(/§ 6. 5. This is known already for the classical symmetry reduction of 
ODEs, which do not respect the reductive algebra structure of the symmetry algebra \Lyij . 

0.3. Discussion. 

Let us make distinction between compatible differential constraints, solvable differential con- 
straints and just differential constraints. The former are based on some explicit compatibility 
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criteria which, provided the type of symbohc syzygy (or type of singularity for characteristic 
variety) is fixeqj, can be written explicitly, see the next section. 

On the contrary there do not exist general solvability criteria for overdetermined systems. All 
of them are based on the prolongation-projection method, which consists in adding compatibility 
conditions and removing irregular points. Thus it is neither explicit nor simple. Indeed the compu- 
tational cost of this process can be quite largelfl. Thus there are no big difference between solvable 
differential constraints and just differential constraints: One starts with an arbitrary overdeter- 
mination and checks consistency getting at the end formally integrable system or contradiction 
(Cartan's formulation for inconsistent equation [CaH as Cartan-Kuranishi theorem states JlaJ. 

In this latter form the method of differential constraints is equivalent to the method of ansatzes. 
It is a matter of taste what to choose then. When we join an additional PDE, it is usually easier 
to check compatibility conditions. On the other hand, an ansatz is ready to provide the form of 
solution, while with differential constraint one should first solve the auxiliary equations. 

1. Details of the method 

Since compatible differential constraints play the central role in reduction of differential equa- 
tions it is important to have efhcient compatibility criteria. The method of E.Cartan, Jane and 
Riquier as well as Spencer theory (see pLV, KLs^ and references therein) provide theoretical base 
and algorithm for establishing compatibility conditions, but they give no formulas. The latter 
were established for a wide class of overdetermined systems of PDEs in |KLi[ |KL2[ IKL4I . 

We briefly recall the results. If £ = {Fi = 0, . . . , = 0} is a system of PDEs (r = dim_ff*^^(£) 
is an invariant of the system expressed via Spencer 5-cohomology) on m dependent functions 
u^, . . . , u™ of n independent variables x^, . . . , x", then we call such system a generalized complete 
intersection iim<r<n-\-m and the characteristic variety and characteristic sheaf over it satisfy 
certain genericity (transversality) conditions (look for details in loc.sit.). The strict inequality 
means that the system is overdetermined, but not too overdetermined. 

For the coordinate situation, which we treat (see IKL4I for a more general context), the multi- 
bracket of (non-linear) differential operators is defined as follows. Let i{G) — {£i{G), . . . ,£m{G)) 
be the linearization of an operator G written in components. Then 

{Fi, . . .,Fm+l} = ^ (-1)" (-1)'^4(1)(-F>(1)) ° • • • 0-^a(m)(-F>(m)) • 

For m — 1 (scalar differential operators) this becomes the Jacobi bracket { , } on diff(l, 1), which 
in the linear case become the standard commutator. 

With the condition of generalized complete intersection a system 8 is compatible (formally 
integrable) iff the multi-brackets vanish due to the system, meaning 

{F,,, . . . , J = mod Jfc,^+...+fc,^^^_i(Fi, ...,Fr) (1) 

for any (m + l)-tuple 1 < ii < • • • < im+i < fi where — ovd{Fj) are orders of the operators 
and Js{F) — (Di o Fi | ordDi + fc^ < s, 1 < i < r) C diffs(m • 1, 1) is the submodule generated by 
Fi, . . . , Fr and their total derivatives up to order s. 

Another important feature is the description of the solutions space. In particular, functional 
dimension and rank of the space Sol(£) can be calculated explicitly IKL3I , which helps to control 
the reduction process, which we shall now discuss. 

^In most geometric approaches to differential equations this is a crucial requirement of regularity, so that the 
type (symbols, characteristics. Spencer (5-cohomology etc) does not change along the equation considered as a 
submanifold in jets. Equivalently one can require that the set of singular points is meager and can be removed. For 
analytic D-modules it can be proved that the set of these point has Zarisky-open complement and thus is negligible. 

^Even polynomial Grcebner basis has double exponential complexity, its differential analogs (,Hu^ and references 
therein) exhibit no decrease, see e.g. IKrI for this question in relation with generalized Lagrange-Charpit method. 
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1.1. Symmetries and Lagrange-Charpit method. 

Let F = {Fi, . . . , Fm)^ be a determined differential operator (this is not only condition on the 
number of functions, but also the claim that determinant of the symbolic matrix does not vanish 
identically). Another operator G = (Gi, . . . , G„i)'^ is a symmetry of F if 

{F,G}^XoF (2) 

for some matrix differential operator A (or its linearization in non-linear casej^. This condition 
means that the overdetermined system F[u] — 0,ut — G[u] is compatible {t is an additional 
independent variable), which can be interpreted as a vector field 9g on the solution space Sol(_F') 
[KLV| (flow of the evolutionary equation) . 

With certain non-degeneracy conditions for the symbols of F, G equation ^ implies that the 
joint system = 0, G = is compatible, which can be interpreted as existence of auto-model 
solutions for G (fixed points of the flow Scjl- Indeed, if m = 1 this was noted already in |KLi| : 
the compatibility condition, specification of ([1]), given by 

{F, G} = A o F + /I o G ( ord(A) < ord(G), ord(^) < ord(F) ) (3) 

is the direct consequence of ([2]). For m > 1 Theorem 1 of IKL4I is not directly applicable, so either 
a version of it for non-commutative rings should be exploited or the reduction via D-modules to 
the scalar system should be applied. This will be performed elsewhere. 

Now overdetermination of the equation F = with another PDE G = (differential constraint) 
is the subject of most reduction methods (since m = 1 in most of them, we restrict our attention 
mainly to this case now). However without transversality condition (complete intersection of 
|KLi| ) the compatibility condition ([3]) will not work. And in fact the reduction may not work 
either, see ^0.11 Thus it is very important to control this general position property (which ad hoc 
is present in all reduction methods, though this is not specified explicitly). 

Classical Lagrange-Charpit method designed to integrate scalar 1st order PDEs consists of 
overdetermining them in a special way. If Fi = is a given scalar PDE, the auxiliary PDEs 
are chosen F2 = 0, . . . , F„ = 0, so that the resulting system is of Frobenius type (i.e. {u^i = 
Fi{x,u)}, I < i < n) and compatible. Two points at this method are to be underlined. 

First is that the condition formulated by them is given in terms of vanishing of Mayer (or 
Lagrange) brackets due to the system, so that the criterion of IKL4I is the higher vector-version 
of the classical theorem. The second is that the symbol of the joint system is of general position 
and the system is not too overdetermined. 

With the general compatibility criterion of fKI^ we can skip the Frobenius condition, but we 
should keep the condition for symbol and that the system is not too overdetermined. The latter 
means that we try to put a minimal number of differential constraints. Of course, this limits the 
general idea, but then an efficient formula for checking compatibility exist^. 

This method was called the generalized Lagrange-Charpit method in IKL2I (it works for vector 
systems as well, with multi-brackets used instead of the Jacobi bracket). The equation G is called 
then an auxiliary integral (the terms differential constraint and side condition can be applied, but 
let us stress importance of non-degeneracy condition). 

Restricting for simplicity to scalar PDEs (this is what basically considered in all methods), we 
can state: 

Most reduction methods are partial cases of the generalized Lagrange-Charpit method. 

This is not true for all methods, since certain ansatzes lead to more overdetermined systems 
and can be interpreted as solvable differential constraints, not compatible differential constraints. 
But nevertheless this makes a convenient specification of the general Olver's claim from 

'^Note that ord(A) < ord(G) for m > 1 and it is < ord(G) for m = 1. 

*At this point I am grateful to S.Igonin and A.Vcrbovctsky for an enlightening discussion of the results of IKL4I 
and the symmetry condition. 

^This is not necessary: Provided that the form of symbolic syzygy (type of singularity of the characteristics) is 
fixed, the compatibility formulas can be specified. This is the essence of the method (not theorem) of [KL4]. 
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1.2. Similarity and other reductions. 

Here we specify how to get major reduction methods from the one considered in the previous 
section. Some of the relations described below are known (see [01] and the references therein), 
but our form is most general and is based on the rigorous compatibility results. We restrict to the 
case n — 2 (independent variables x,t), m — 1 (dependent variable u) and one single PDE (most 
reduction methods are discussed precisely for this particular case), though versions for the general 
case are straightforward. 

1. Symmetries and generalized symmetries. The point symmetries are quasi-linear auxiliary 
integrals of the form: 

G — a{t, X, u)ut + b{t, X, u)ut + c(i, x, u). (4) 

The only difference between classical symmetries introduced by S.Lie [Liei ILE| and generalized 
symmetries of Bhiman and Cole |BC| , see also |0v| IMej , is that they satisfy different relations ([2]) 
and ([3]) respectively, of which the latter is more general. Thus for establishing (partially) invariant 
solutions one can treat them without any distinction. 

Contact symmetries were, also introduced by S.Lie and they are general 1st order constraints 
G = G(t, x,u,Ut,Ux), which satisfy ([2]). Though never treated generalized contact symmetries 
should be the general first order G satisfying ([3]). They are more plentiful than generalized point 
symmetries. 

2. Higher symmetries. These are just differential operators G of higher order satisfying 
([2]). Upgrading it to ([3]) we get higher generalized symmetries or auxiliary integrals (as well 
as compatible differential constraints and side conditions). The only restriction which should not 
be forgotten in all these overdeterminations is the non-degeneracy condition (i.e. _F = G = is a 
complete intersection). 

Higher order auxiliary integrals exist in abundance, but to find from their continuum certain 
concrete families one should specify the form of G due to original PDE (see examples in |KLi| ). 

3. Direct method and beyond. Direct reduction of Clarkson and Kruskal |CKj is an ansatz for 
the form ([4]), so it gives the same exact solutions as explained in 1. This method was generalized 
by Galaktionov to the following ansatz: u — U{t, x, wi {({t, x)), . . . , WkiCit, x)) (he took k — 2) and 
the substitution is aimed to reduce the PDE to an ODE of higher order for the variable z = C,{t, x). 

If V is the vector field with first integral C{t,x), then Olver [Ol] demonstrates that this is 
equivalent to a differential constraint of the form G — V'^{u) — $(t, x, u, V{u), . . . , V'^^^{u)) — 0, 
which is a special form (degenerate symbol, characteristic variety is fc- multiple point in <CP^) of 
quasi-linear PDE. 

Provided that differential operator F of order m has Char'''(F) not containing this point and 
that {F, G} = mod Jfc+™_i(F, G), the PDE F is reduced to an ODE of order km (KL^]. 

Remark 3. If an ansatz reduces PDE F = of order k to an ODE of order N or a system with 
solutions space of dimension N , and number N is not divisible by k, then this reduction cannot 
be obtained by a single compatible differential constraint. Indeed transversality is necessary for 
reduction and then compatibility implies dimSol(i^, G) = ord-F-ordG. Thus this case corresponds 
to several (different orders) jointly compatible differential constraints or to solvable differential 
constraints (this can be a single scalar G). 

Generalization of this method to arbitrary n and m is this: One specifies G to have degenerate 
form and characteristic variety. Darboux integrability then becomes a partial case. 

4. Separation ansatzes. Additive separation ansatz u ~ f{t) + g{x) (we continue the case 
n = 2, m = 1 for simplicity) correspond to G = u^y^ multiplicative one u = f(t)g{x) corresponds to 
G — uutx — UtUx- Functional separation u — U{f{t)-\-g{x)) generalizes these both and corresponds 

to G = Utuluttx - ujUxUtxx + UtxiulUxx " ulutt) ■ 

Functional separation of variables u = X^iLi fi{^)9i{^) corresponds to a complicated non-linear 
higher order constraints. Note that for fc > 1 encoding the ansatz cannot be chosen as a single 
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PDE. For instance, when fc = 2 we can have either a mixed ansatz-constraint (i.e. an auxiUary 
integral G containing an unknown function) 

{utxU - UtU^)p"{t) + {u^Utt - UuxU)p'{t) + {uttxUt ~ UttUtx)p{t) = 

or 2 compatible PDEs of order 5, expressing non-trivial solvability of the above equation for p{t). 

5. N-soliton and finite-gap solutions of soliton equations. These are always special type 
constraints. Consider for instance KdV equation ut — u^xx + QuUx- iV-soliton solution is obtained 
from an ansat^I^ u = cwxx, where e"' = ^ Cie"'*+'''^. Functions w of this type can be described 
as solution to a PDE F[w] — e^™i[e"'] — 0, where i is a linear differential operator with constant 
coefficients, is a non-linear operator and there exist other non-linear operators of high order 

such that ^ o F[w] = ^ o cV^. Then the (solvable) differential constraint is <i>[u] = 0. 

The space of A^-gap solutions (with specified place for gaps) is a finite-dimensional (symplectic) 
manifold invariant under KdV dynamics. The standard description is however in non-differential 
terms (spectral properties of the Hill equation) , but this can be re- written via differential constraint 
with r-functions as coefficients. Since the latter are in turn solutions to differential equations, the 
auxiliary integrals can be specified in an elementary, but very non-linear form. 

Notice also that Painleve test is a special kind of ansatz and so can be seen as a constraint. 

6. Evolutionary equations and inertial manifolds. The generalized Lagrange-Charpit method 
is much simpler for evolutionary PDEs. Indeed li F — ut — $(t, x, u, Ux, Uxx, ■ ■ so that no t- 
derivative enters into the last term [x and u can be multidimensional), then an auxiliary integral, 
with no loose of generality, can be represented in the form G = ^'(t, a;, u, Ux, Uxx, . . .). Indeed all 
derivatives of ut can be expressed from the equation F = Q. 

Seeming difficulty that the characteristic variety of the system f = {F = G = 0}is degenerate 
for 1-dimensional space variable x and ord<i> > 1, can be resolved by introducing higher weight 
for the variable t |KLi| . If we assume moreover that the auxiliary integral does not contain t 
explicitly, i.e. G = G{x, u, Ux, Uxx, ■ ■ ■)^ then G = gives an equation of the stratified submanifold 
iS (if n — 1 = dim{(a;)} > 1 it is an infinite-dimensional manifold, while in the case n = 2 we get 
an ODE), which is invariant under evolutionary flow, see [KH ILLj for examples. 

Such (finite-dimensional) invariant submanifolds S = {G = 0} are called inertial if in addition 
they attract evolutionary dynamics ut — In this case E contains an attractor, which can be 

investigated by the method of dynamical systems. When exist, inertial submanifolds can be found 
via compatible differential constraints. 

1.3. Examples. 

Since the method is the most general, we cannot classify all solutions obtained via it (this would 
be just the description of Sol(f )), but we show some relations with what other did and something 
new as well (see also |KH ILL) for other examples). 

We will restrict to the case n = 2, m = 1 and one PDE F — Q. Assuming complete intersection 
condition, its compatibility with the differential constraint G = is given by [F,G]s = 0, where the 
latter bracket is the reduction of the Jacobi bracket {F, G} by the subspace Jord F+ov<iG-i{Fi G); 
this reduction is called the Mayer bracket |KLi| . 

1. For F — Alt — 1 (where A = -\- is the Laplace operator on R^) the constraint 
G = (%/«)* + — a is an auxiliary integral iS a — ^ ot j. This follows from the formula 

[F, G]£ = (l-20Kl-40) _ ^Q^g ^^^^ Q 

is non- linear, while F is linear non- homogeneous (cf. i i2.2p . 

2. For F = c{t)utt - Uxx we let G = Ut + a{t)u. Since u ■ [F,G]£ = c{t)a"{t) + c'{t)a'{t) - 
2c{t)a{t)a' [t) ~ a{t)'^c'{t), the compatible constraint of the given form corresponds to the function 
a{t) — Co + ci J -^jj. It is not a point symmetry, so it is a non-classical symmetry. 

3. Let us re-expose an example from [OR . For the heat equation F — ut — Uxx they take 
G — tut — xux + ix^ . This is not a compatible differential constraint, because H — Gjg = 

^"The parameters should be in fact more specific, otherwise this family yields more solutions. 
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Ut — 6x, which is not identically zero. But adding this PDE makes the system compatible. Indeed, 
G = H = is a compatible Frobenius system with the solution u = + 6tx + c G Sol(i^) . 

4. We continue with the heat equation F — ut — u^x- The constraint of the multiplicative 
ansatz G = UtxU — Utu^ does not satisfy the complete intersection condition. The compatibility 
condition is simpler: H = u'Dt{F) + 'Dx{G) — UtF = uuu — = 0. Thus again the system 
F = H = is not compatible, but the "projected" system F = G = H = is oi Frobenius type 
and is compatible. 

5. Finally let us consider a more complicated case. Choose 

F = Ut - Uxx - f{u) = 0. (5) 

This is a non-linear reaction-diffusion equation (called also Kolmogorov-Petrovsky-Piskunov and 
non-stationary heat equation) . The following is an auxiliary integral: 

G = UxUxxx - ul^ - f{u)uxx + f'{u)ul = 0. (6) 

This follows readily from the identity [F, G]£ — 0, but can be computed directly from invariance. 

Now an important thing is that PDE ^ was a popular target for investigation for exact 
solutions. Detailed analysis and a review is contained in [CM] . In all the results f{u) has some 
particular form, most often cubic. In |KL] this equation was investigated with a different method 
and still some exact solutions for some other particular form of /(u) were obtained (an exception 
was with first order auxiliary integrals, when f{u) was general, leading to traveling waves). 

But the auxiliary integral ^ works for all f{u) and moreover it is already an ODE for u in 
X, so it is the reduction (the solution of the evolutionary equation is reduced to two quadratures: 
first is G = and the second is given by integrating the vector field corresponding to F on Sol(G), 
see pCLllLLj V 

Remark 4. This shows that the method of differential constraints (Lagrange- Charpit in this case) 
is more powerful than the ansatz approach. The reason is that substitution of separation ansatzes 
leads to a functional equation, which is difficult to solve in general (indeed following Abel we solve 
functional equation via differential ones). 

6. In order to explain the result of 5 let us mention that ([5]) was investigated for functional 
separation of variables (as in direct reduction method) u — U{(p{t)x + ip{t)) in [PZ| . The corre- 
sponding functional equation was partially solved with certain type of f{u) specified. 

Now the above ansatz is equivalent to the constraint ( — ) = 0, which coupled with ([5]) yields 

^x 

a forth order PDE (after multiplication by u^) 

Gi = uluxxxx - i^UxUxx + f{u)Wx)Uxxx + 2u^2: + '^fi'u)ulx ~ f'{u)ul.Uxx + f"{u)u')^ = 0. 

However it is not compatible with (O: [F, GJg = 2 G • G2, where G is given by (O and G2 = 

f{u)UxUxxx - 2f{u)ul^ + f'{u)ulUxx " f"{u)u'^^. 

So the system splits. The first branch = G = is compatible and was noted above. This 
choice gives unconstrained f{u). 

The other system F = G2 = is not compatible. Compatibility analysis implies that the 
system F = G2 = is solvable iff / = f{u) satisfies the following highly nonlinear ODE: 

2f\f'ff(^) + f^f'jm^^ff + i2//(4)) + //'/(3) (3o(/")4 + 45//"(/(=^))2 + 16/(/")V(') 
+ 12/2(/W)2) + /'2 (3(/")^ + 39/(/")^(/(^))^ + 3/(/")V(') + 18f (./(^))V'*^ + 4/V"(/W)2) 
+/(8(r)'+31/(/'r(/(^))^+45/V"(/^'^)V(')-6/^(/")'(/''^)'+/'(-27(/(3))^-}-2/(/W)3)) ^0 

Thus in this case we also get a reduction, but again for a special kind of f{u). 
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7. Note that the operator in ^ is F = dt — L with L[u] = u^x + f{u) and G = V o L with 
V = UxDx — Uxx- Thus the evolutionary flow of F is given by the second order ODE on the 3- 
dimensional space Sol(G) of solutions of ODE G = in x-variables only. Note that Sol(L) C Sol(G) 
is the 2-dimensional surface, on which the flow is stationary. On the complement the flow is non- 
trivial. This geometric reasoning shows that while there are many auxiliary integrals of the form 
G = DoL, not every such G is a compatible constraints, but any is a solvable differential constraint 
(indeed Sol(L) wiU be a part of Sol(F, G)). 

Let us also note that the solutions of (l5|)-|-([6]) can be obtained by a Frobenius overdetermination. 
Indeed it is equivalent to the joint system 

Utt ^ -k{ut ~ f{u)), Utx = {Ut-f{u)), Uxx=Ut-f[u), 

Ui Ux 

which is compatible (F = is the last equation). 

1.4. Algebraic structure of the set of auxiliary integrals. 

In many papers devoted to differential constraints (for instance |01[ IClj ) it is noted that the 
method is too general, so it is impossible to operate with, contrary to the symmetry group ap- 
proach. It is indeed quite general due to Theorem [TJ But being reasonably restricted it remains 
quite general/powerful, while will possess a certain algebraic structur j^. which we are going to 
explain now. 

For simplicity we restrict to the case of a single scalar PDF F — 0, where F G diff(l, 1) is a 
(non-linear) differential operator over M, as this is an often target of investigations (but change 
of the bundle 1 to tt is possible with some assumptions on the operators) . 

Remark 5. Non-classical symmetries are special kind of auxiliary integrals in generalized Lagrange- 
Charpit method, so their totality does not possess algebraic structure as well. Sometimes to stress 
group approach it is said that generalized symmetry G of F — can be characterized by the prop- 
erty that the flow of it preserves the solutions Sol(i^, G) (contrary to the larger space Sol(i^) for 
classical symmetries). But then general auxiliary integrals are the same kind of generalizations of 
the higher symmetries (we can treat the flow as defined formally). 

In what follows we write £f = '^Fpa'^a the operator of linearization. Here {x^,pc) are the 
canonical coordinates on the jet-space J°°{M) induced by a coordinate chart on the base M", 
a = («i, . . . ,z„) is the multi- index and Va- — 'D\^ ■ ■ • 2?^" is the higher total derivative operator 
with TJi = + ^Pa+udp^ ■ Thus the operator of universal linearization transforms non-linear 
operators diff(l, 1) to linear operators Diff (1, 1) if we fix a jet Xoo G J°^{'^, !)■ More generally it 
is the map 

i : diff(l, 1) ^ Diff(l, 1) = G°°(J°°(1)) ®c~(A/) Diff(l, 1) 

(the latter are called scalar "^-differential operators). 

Let us denote Sym(F) = {H e diff(l, 1) : {F, H} = £gF for some 9 e diff (1, 1)} the space of 
symmetries of F, which is a Lie algebra w.r.t. the Jacobi bracket. The defining relation determines 
the family of vector subspaces Symg{F) = {ff G diff (1, 1) : {F, H} = £gF}. Note that Symo(F) is 
the Lie subalgebra of Sym(F) and Symg(F) is the module over it (considered as a vector space!). 

Notice that the following space is a vector space, but in general it is not a Lie algebra: 

sym*(F) = {9} = Sym(F)/ Symo(F). 
Next for /i e diff(l, 1) with ord(/i) < ord(F) denote 

Aux^(F) = {G e diff(l, 1) : {F, G} = £xF -f £^G for some A G diff(l, 1)} 

^^In too general form, we observe no algebraic structure for compatible differential constraints. The same concerns 
solvable constraints, even if the form is prescribed. 
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and let Aux(_F) = U^.ord(/^)<ord(_F) Auxp(i^) be the set of all auxiliary integrals. Notice that we 
omitted here the condition of general position for the symbols of F and G (complete intersection) . 
With this condition the joint system {F = G = 0} is compatible and G is a compatible differential 
constraint. 

Aux^(F) is a vector space. Changing G i-^ G + trjF we change coefficient A in {F^ G}. The 
equivalence classes form the set aux*(i^) = {/i}, which in general is just a topological space. 

It's easy to see that Aux(i^) is neither a vector space nor a module over Sym(F) (not even 
over Symo(i^) as vector space). Indeed, if G € Aux(F), {F, G} = ixF + ^^G, and H G Symg(F), 
then {F, G + H} ^ ix+gF + ^^G. Only Auxo(F) = Sym(V) is the module. 

However it has the structure of algebroid type. To describe it we need to introduce some 
notations. Denote ad^ = {H,-}. Let 9c be the evolutionary differentiation, which is the dual 
operator to universal linearization |KLV] : 9gF = ipG. Note that adfi — £h — 9h- 

The Hessian operator Hess : difr(l,l) x diff(l,l) '^^Diff(l,l) is defined by the formula 
Hessi? G = [9g,^f]- In canonical coordinates we have: 

HessF G{H) - HessF(G, H) =Y, Fp^pMG)VriH). 

For linear differential operators F the operator Hess_F vanishes. These operators satisfy the fol- 
lowing properties (direct check): 

Proposition 2. 1. Symmetry: HesspiG, H) = HesSi?(iJ, G). 

2. Compensated Leibniz rule: {F, ioH} = ^{f,g}H + iclF, H} - Hess_F(G, H). 

3. Linearization anomaly: [tpi ^g\ — ^{f,g} + (Hessf? F — Hess^? G). 

Now we can use this to calculate: 

{F, {G, H}} = {{F, G}, H}} + {G, {F, H}} = -{H, l^F + l^G}} + {G, leF} 
= {^{XM} + [txJe] - he.G} + HessH A - Hessg e)F + £^{G, H} 

+ (^{p,//} - + HessHA*)G = -{adH+£ey^G modF. 

From the last two lines we conclude two actions: 

Proposition 3. 1. If H e Symg(J^), (^ad„ -f^e - Hess^ )Ai = 0, then adj? : Aux^(F) Aux^(F). 
2. IfHe SyiJig{F), i.e {adn +te)F = 0, and [adn +ee)e^ = 0, then adn : Aux^(i^) ^ Sym(i^). 

Thus we get partial algebraic structure on the set Aux(F), so that the fibers of the vector 
bundle Sym(i^) sym*(F) act on the fibers of the vector bundle Aux(F) aux*{F). 

2. Other related results 

The compatibility criterion via multi-brackets discussed in the previous section can be gener- 
alized to the case, when we still have r < m + n equations on m dependent functions, but allow 
degenerations keeping regularity. Let us begin for simplicity with the case m = l,r = 2. For 
complete intersection we required that the characteristic varieties Char'''(i^) and Char*'(G) are 
transversal in each fiber CP""^ ( |KLi| ). We can relax this by letting them have common compo- 
nents. Assume the symbols a(F), a(G) have a (greatest) common divisor q. Then we can reduce 
the bracket by an operator Q with a{Q) — q. 

In fact, let k = ord(F), I = ord(G), t — ord((5) and S, T be differential operators with symbols 
a{F)/q, a{G)/q (regularity assumption means the division is made over all points of the system 
£ — {F — Q,G — {)} and that loci of polynomials a{F)/q,a{G)/q intersect transversally) . Then 
£ is formally integrable iff £tF — isG — modJ^k+i-q-i{F,G). The operator on the left is the 
reduced analog of the bracket. 

Generalization to the case 2 < r < n is straightforward. Multi-bracket analog m > 1 is 
more involved. It is also possible to have explicit criteria when r > m + n, but this will be treated 
elsewhere. In this section we consider some questions related to compatibility issues of establishing 
exact solutions. 
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2.1. ODEs. 

In [Cl]§6.3 it was noticed that contrary to classical symmetry reduction the nonclassical and 
direct method are difficult to use for ODEs. In fact, there's no difficulties if we view them as 
differential constraints. However in this case the number of independent variables cannot be 
reduced (it is already 1) and the order should be reduced instead (cf. Remark [5]). 

Consider an example, which has a general nature. Let 

u" — F(x, u, u) (7) 

be an ODE. Seek for the compatible differential constraint in the form u' ~ if(x^u). In this case 
the second equation is just the intermediate integral of the first equation and condition of this is: 

(fix = ~Vu'P + F{x,u,ip). (8) 

Thus reductions of ODE ([7]) are bijective to partial solutions of PDE This is an analog of 
Hamilton- Jacobi approach. 

It is not clear how to find a solution to PDE, but whenever it is done we proceed to simplify 
([7|). For instance all classical reductions can be easily recognized. Let's mention the basic: 

1) If F is independent of x, then x does not enter ([5]) and we seek for solution = ip{u). Then 
© is the ODE-reduction. 

2) If F is linear in u, then from ([8]) we read that it worth searching Lp to be homogenous in u 
of degree 1. 

3) If F{e°'*x, el^^u, e^'^'^^V) = e^'^-^°''>^F{x, u, ip), one can take ip{e°'*x, e'^'u) = e^'^-°''>*ip{x, u). 
In all these cases one uses a symmetry of ([7]) to be the symmetry for ^ , which in the considered 

cases are respectively: dx, udu and axdx + (3udu- 

On the other hand what we do is just to solve ([8]) with the help of symmetries, which we can 
write as auxiliary equations: 

1) ipx = 0, 2) uipu - </3 = 0, 3) aipx + /3(^u + (a - f3)<fi = 0. 

But in general to solve ([8|) we can use more general Lagrange-Charpit method with some auxiliary 
equations tpu = G(x, u, ip, ipx), so that the procedure of solving ^ becomes repetitive. 

2.2. Linear systems. 

Linear PDEs, as the simplest equations, produce sometimes confusion in problems related to 
compatibility. Here we discuss some. 

1. Olver claims in Theorem 13 of [Olj that a compatible differential constraint for a second 
order linear equation must be linear. This is wrong. 

A simple counter-example can be obtained for the most popular equation of the sort: heat 
equation F = Ut — Uxx- Then G — {ul + u^) cosh(2x) — 2utUx sinh(2x) -I- 2e^*(itt — u) is an auxiliary 
integral, i.e. the system £ — {F = G = 0} is compatible. 

The order of G is 1 and Sol(£) = {cie*^^ -fC2e*+^ +c\ -\-c^}. The point is that this is a curved 
surface in the linear space Sol(i^), in fact a quadric. By this reason the function G is quadratic. 

Similarly if we take a 4-dimensional curved surface in Sol(-F), we get an auxiliary non-linear 
integral G, i.e. a compatible constraint contradicting Olver's claim (he was concerned with second 
order nonlinear constraints) . For instance the following auxiliary integral works: 

G — a;^(2u( -f 2U(( — uu + ut) + 2utt(l — t -\- 2xutx — 2xux) 

+ 2ut(t ~ 2x^utt - 2xutx + 2xux) + 2utx{utx - x - 2ux) + 2(xux + Ux"^ - u). 

It is quite clear that the construction is most general and plenty of nonlinear constraints for 
linear equations exislF^. For instance F — uu^Uxx possesses the compatible differential constraint 
G = detHess(u) -I- ae''" with Monge- Ampere type non-linearity (see |KLi| ). 

Constant coefficient linear PDEs are often target of investigation for particular solutions via reduction methods 
to excitement of applied mathematicians, who find all solutions by the Fourier transform. 
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2. For linear PDEs with constant coefRcients an auxiliary integral can be chosen of the same 
form: A system of PDEs with constant coefficients is compatible, provided it is of generalized 
complete intersection type, as follows from the compatibility criterion of Section [TJ Without this 
requirement the claim is wrong. 

For instance the system {utx = u,ut = u\ has only trivial solution (m = is always a solution of 
a linear homogeneous system) . Also for linear non- homogeneous system with constant coefficients 
the above claim would be wrong: The system {ut = u -\- \,Ux — 1} has no solutions. 

Since any two linear operators Li,L2 with constant coefficients commute {Li,L2} — 0, one 
is a symmetry for the other. But without complete intersection condition the PDEs Li[u\ = 0, 
L2[u] — might be incompatible. This shows importance of our generic position assumption. 

3. Investigating whether nonclassical method gives more symmetric solutions than the classical 
one, when applied to linear PDEs, the authors of j AGBj make the following statement: 

Every solution of a linear PDE is invariant under some classical symmetry. 

Then they prove it under the following assumptions: a) n = 2 (two independent); b) m = 1 
(one dependent); c) the PDE has constant coefficients. 

It is a bit disappointing that the proof contains determining relations for symmetries, calcu- 
lations etc. As the authors expect that their theorem can be extended to bigger number n, we 
confirm this and make the following generalization: We will consider systems of r equations on m 
unknowns and will not require that coefficients of the PDEs be constant. 

Proposition 4. Let L G Diff(TO • 1, r • 1) be a linear operator. Then every solution to the system 
L[u\ = is invariant under some symmetry provided the following condition on the symmetry 
algebra holds: 

Sym(i) ^Sym,(i)~Sol(i). 

Here we denote by Sym^(i) C Symo(i) the subalgebra of shifts by solutions. The symmetries 
have the form gdu ■— ^g^{x)duo for g = {g^,...,g"^) e Sol(L). Thus such a field has the 
generating function lifted from the base |KLV| and its prolongation is ^ Da{g^)d^] . The inequality 
above holds true for operators L with constant coefficients {d^i is a symmetry, so that we can 
allow dependence on x^, . . . , x") and thus the claim contains Theorem 1 of [AGBj . 

Proof. Let w be a solution and ^ a symmetry from Sym(L) \ Sym^(L). Denote g — S^u — u. 
Then g S Sol(L) and the symmetry ^ — g9„ leaves u invariant. □ 

The large size of symmetry algebra does not help in solving linear PDE because the piece 
Sym^(L) is generally untractable. One should rather consider the finite piece Sym(L)/ Sym^(L) 
(it is given by linear overdetermined system and often is finite dimensional). 

Thus even though the conclusion of [AGBj is that no generalized symmetry gives new solutions 
compared to classical symmetries, this claim does not bear a constructive obstruction for consid- 
ering generalized symmetries (compatible constraints) . For me though footnote [T^] does, when 
PDEs with constant coefficients are considered. 

2.3. Conclusion. 

In this paper we demonstrated how one can establish new exact solutions of non-linear PDEs. 
The scheme we describe is general enough to include most practiced methods. The method is 
indeed classical, but we furnished it with an efficient underlying machinery of investigation for 
compatibility and precise general position conditions to hold. 

The generalized Lagrange-Charpit method we discuss was essentially developed in IKL2I and 
here we make a modification, allowing non generalized complete intersections, but with a fixed 
type of the syzygy module (then a criterion different from that of |KLi| may show up) , or changing 
compatibility claim to solvability, of which vanishing of multi-brackets due to the system IKL4I is 
the first and most complicated step. 
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Thus the proposed technique is not too general, but manageable. We have not performed 
too many calculations in this paper (see |KLi[ IKL[ ILL] ), but frankly explained how to do it, 
leaving more results for future publications. In addition, the method allows to satisfy physical 
pre-requirements for the desired solutions (form of the compatibility constraint), which might be 
too restrictive with symmetry approaches. 
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